Weekly <- readRDS("Alergia.rds") %>% group_by(Semana) %>%
summarize(mean = mean(platano_oriental, na.rm = TRUE),
max = max(platano_oriental, na.rm = TRUEl),
min = min(platano_oriental, na.rm = TRUE),
sd = sd(platano_oriental, na.rm = TRUE))
ggplot(Weekly, aes(x = Semana, y = mean)) + geom_ribbon(aes(ymax = max,
ymin = min, fill = "red")) + geom_ribbon(aes(ymax = mean +
sd, ymin = mean - sd, fill = "blue"), alpha = 1) +
geom_line() + scale_fill_manual(name = "leyenda",
values = c("blue", "red"), labels = c("Error estándar",
"Extremos")) + ylab("polen de platano oriental /m³ de aire") +
theme_classic() + theme(legend.position = "bottom") +
scale_x_continuous(breaks = seq(from = 2.5, to = 49.5,
by = 4), labels = c("Ene", "Feb", "Mar", "Abr",
"May", "Jun", "Jul", "Ago", "Sep", "Oct", "Nov",
"Dic"))